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Abstract 

We study the dynamics of a phantom scalar field dark energy interacting with dark matter 
in loop quantum cosmology (LQC). Two kinds of coupling of the form ap m (f> (case I) and 
3/3H(p ( f > + p m ) (case II) between the phantom energy and dark matter are examined with the 
potential for the phantom field taken to be exponential. For both kinds of interactions, we find 
that the future singularity appearing in the standard FRW cosmology can be avoided by loop 
quantum gravity effects. In case II, if the phantom field is initially rolling down the potential, 
the loop quantum effect has no influence on the cosmic late time evolution and the universe will 
accelerate forever with a constant energy ratio between the dark energy and dark matter. 
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I. INTRODUCTION 



Recently, present accelerating expansion of our universe has been confirmed by many 
observations, such as the cosmic microwave background (CMB) anisotropy, type la super- 
novae, and large scale galaxy surveys [H, [2| • In order to explain this observed phenomena, 
dark energy is assumed to exist in the universe within the framework of general relativ- 
ity. Dark energy is an exotic energy component with negative pressure and accounts for 
about 72% of present total cosmic energy. In addition, observations also show that there 
is another dark component in the universe, i.e., dark matter, accounting for about 25% 
of total cosmic energy today. The simplest candidate of dark energy is the cosmological 
constant. It is, however, plagued with the so-called coincidence problem and the cosmo- 

i r 

logical constant problem [3]. Thus some dynamical scalar fields, such as quintessence [4j, 
phantom quintom [f| and hessence Q, are proposed as possible candidates of dark 
energy. It is worthy to note, however, that for these scalar field models the coincidence 
problem still remains. Although the two dark components are usually studied under 
the assumption that there is no interaction between them, one can not exclude such a 
possibility. In fact, researches show that a presumed interaction may help alleviate the 
coincidence problem [s| . Therefore interacting dark energy models have attracted a great 
deal of interest. 

Most of the present studies on dark energy are carried out in the framework of classical 
Einstein gravity. However, it is commonly believed that quantum gravity effects would 
play a role in the evolution of the universe. Therefore, it is desirable to examine the 
properties of dark energy in a theory of quantum gravity. One of such theories which we 



are interested in the present paper is Loop Quantum Gravity (LQG) (see e.g. P. [lol.llll.ll2| 
for reviews), which is a nonperturbative background independent theory. At the quantum 
level, the classical spacetime continuum is replaced by a discrete quantum geometry and 
the operators corresponding to geometrical quantities have discrete eigenvalues. LQG has 
been applied in cosmological context as seen in various literature where it is known as 



Loop Quantum Cosmology (LQC) [13, LL4 ll5|, ll6|]. The effects of loop quantum gravity 



modify the standard Friedmann equation by adding into it a correction term — p 2 / p c which 
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essentially encodes the discrete quantum geometric nature of spacetime 

QHQ- Let 

us note that although extra terms such as those derived in Ref. 19j are in principle 
possible, the effect of such terms is negligible for all practical purposes 20]. When this 
correction term becomes dominant, the universe begins to bounce and then oscillates 
forever. Therefore both the future singularity and the singularity at semi-classical regime 



can be avoided 
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,122, 



231 ] . Recently, the dynamics of phantom, quintom and 



23|, 
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251 ] and behaviors different 



hessence in loop quantum cosmology have been studied 
from that in the standard FRW cosmology, such the avoidance of future singularities, are 
found. At this point, it should also be noted here that the LQC modification, although 
very interesting, is only derived for pure isotropy and is not very clear how it fits within 
broader models. 

In this paper, we will discuss the dynamics of a phantom scalar field dark energy, 
coupled to dark matter in loop quantum cosmology, to see in what way the LQG effects 
would affect the cosmological evolution of the system. Two kinds of interactions of the 
form, acp m <j) 26] and 3(3H(p ( j ) + p m ) jsj], will be studied. It is worth pointing out that the 
first coupling could arise in string theory or after a conformal transformation of Brans- 
Dicke theory, while the second is motivated by analogy with dissipation of cosmological 
fluids and has been proposed for a possible dynamical solution to the coincidence problem. 
Let us note here that the phantom field, although a viable candidate of dark energy, also 
has some strange properties, such as being unstable to vacuum decay and the violation 
of the dominant energy condition. 



II. LOOP QUANTUM COSMOLOGY 



By incorporating the effects of loop quantum gravity which essentially encode the dis- 
crete quantum geometric nature of spacetime, the effective modified Friedmann equation 
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in a flat universe is given by 



3, 0, 0, 
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p 



Pc 



p 



(1) 



where H = a/a is the Hubble parameter, p is the total energy density, and a dot denotes 
the derivative with respect to cosmic time t. Here we set 8nG = 1; and the critical loop 
quantum density is 

V3 



where 7 is the dimensionless Barbero-Immirzi parameter. Let us note here that it has 



been suggested that 7 ~ 0.2375 by the black hole thermodynamics in LQG) [21 
Differentiating Eq. ([1]) and using the conservation equation of cosmic total energy 

p + 3H (p+p) = 0, 

one obtains the effective modified Raychaudhuri equation 
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H 



1 



-(p + p)(l-2^ 



(3) 



(4) 



where p is the total pressure. Actually, as shown in [29( , the effective modified Raychaud- 
huri equation can be also derived directly by using the Hamilton's equations in LQC, 
without assuming the energy conservation. 



III. DYNAMICS OF THE INTERACTING PHANTOM SCALAR FIELD DARK 
ENERGY IN LQC 

Let us suppose that there are only the phantom scalar field dark energy and dark 
matter in a spatially-flat universe. The Lagrangian for the phantom scalar field fluid is 

£ = (1/2)9^0-^(0), (5) 

where V(4>) is the potential of the phantom field. Therefore, the energy density and 
pressure for the phantom field can be expressed as 

P<t> = ~\4> 2 + V(0) , (6) 

1 It is interesting to note that this kind of modified Friedmann equation also appears in cosmological 
braneworld models with a single timelike extra dimension [28[ 
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and 



]<P 2 - V(4>) . (7) 



We assume that there is an interaction T between the phantom dark energy and dark 
matter. A positive T corresponds to energy transferring from phantom to dark matter 
and vice versa for a negative one. Therefore the energy densities for the phantom scalar 
field and dark matter satisfy the following equations 

p^ + 3H(l + w^)p^ = -r, (8) 
Pm + 3H Pm = +r, (9) 

and the Raychaudhuri equation. The evolutionary equation for the phantom field and the 
modified Fiedmann equation can be expressed as 

/r = -i(, m -0*)(i-^), (io) 

<j> = -3H4> + V' + ^, (11) 

H> = !(«.- ^ + V) (l - £) , (12) 

where V = dV/d<p. In this paper we only consider the case of exponential potential 
V(4>) = Voexp(— A0) with a positive constant A. 

In order to study the dynamics of the above system, we introduce the following dimen- 
sionless variables 



_ _ VV = P_ ±_ = ]_d_ 

X ~V6H J V ~ V3H ' Z ~ Pc ' dN~ Hdt' 1 J 

where N = In a is the e- folding number and is used as an independent variable instead of 

cosmological time. Using the above definitions, the effective modified Friedmann equation, 

namely Eq. (Tl2|) . can be rewritten as 



( Pi 



and Eq. (TlOi) becomes 



H 



\3H 2 



3 



x 2 + y 2 (l-*)-l = 0, (14) 



2 V 1 



+ x 2 - y 2 ) - 3s 2 



[l-2z). (15) 
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In addition, we will be interested in three scalar quantities, which are, respectively, the 
phantom fractional density parameter fi^, the effective phantom equation of state and 
the effective equation of state for total cosmic energy w e ff given by 



fL, = -x 2 + y 2 



x 2 + y 2 

9 9 : 

x z — y z 



2 2 

x -y 



Using the Eqs. (fTT| [TH [T5|) . we obtain the autonomous equations: 



'3 r 
x' = —3x — \l-Xy 2 -\ — 



— x 



3x< 



1 



2 Xx v - y 



3x z 



2 \1 - z 
-3z-3z(l-z)(-x 2 -y 2 ). 



+ x 2 - y 2 



2 VI- z 

1-22) 



i 2 2 

+ x — y 



(16) 



(1-22), (17a) 

(17b) 
(17c) 



Let / = x', g = y',h = z'. Then using the following condition 



(/, 9, h)\ 



(x c , y c , z c ) 







we can obtain the critical points (x c , y c , z c ) of the autonomous system. 

Next, we will study, by examining the stability of the these critical points using the 
standard linearization and stability analysis, the dynamics of two different interaction 
cases, i.e., I\ = ap m (f) and T 2 = 3j3H(p m + p^) between the dark energy and dark matter. 



A. Case I: I ap m (j) 

This kind of interaction could arise from string theory or scalar-tensor theory [26]. 
The dynamics of this interacting phantom model was studied in the standard FRW cos- 
mological framework in Ref. [30( and it was found that energy transfer either from the 
phantom field to the dark matter or vice versa yields similar cosmological consequences, 
and the energy density of the phantom field increases with the cosmic expansion, leading 
to unwanted future singularity. In this section, we will study the cosmological evolution 
of this interacting model in the framework of loop quantum cosmology. The autonomous 
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Eq. (flTI) can be rewritten as follows 



x 



• ■ 2 / 1 2 2 

-6x — \ -\y —a I h x —y 



—x 



3x 
3 



2 V 1 - z 



1-z 



, 2 2 

+ x — y 



(l-2z) 



Axy - y 



1 



2 VI - -2 

2 „.2 



i 2 2 

+ x — y 



(1 - 2*) 



-32-3z(l -z)(-x 2 -y 2 ) 



(19a) 

(19b) 
(19c) 



Then we obtain five critical points: 



Point (A) : {^-, 0, 0) 



Point (B) : 
Point (C) : 



A 2 



A 

71' 

A / A^ 



FVi„l (D): 0, 

za 3 



, v 7 ^ 1 V2a 2 - 3 + 2aA 
Point E : - -, , 



(20) 
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The eigenvalues, /x, of the coefficient matrix of the linearized equations for these critical 
points A, B, C, D and E can be expressed respectively as 



Point (A) : 

" 2 



, 2 3\ 2 2a 2 -3 + 2aA 
Hi = - [a + - , /i 2 = 2a -3, /i 3 = . 



Hi = A 2 , /i 2 = -3-iA 2 , /i 3 = -3 - X(a + A), 



Point (B) : 
i = A 2 , 

• Point (C) : 

/ii = A 2 , /j 2 = — 3 - ^A 2 , // 3 = -3 - A(a + A), 

• Point (D) : 



3 / / 6 - 4a 2 



3 / /\ 6 - 4a 2 

3A 
2 a 
•Point (E) : 

-3A 

Mi 



a + A ' 



3 2a + A/ / 8[3 + A(A + a)][2a 2 + 2aA-3] \ 
^""I^TX^ + V 3[2a + A] 2 J' 



3 2a + A/ / 8[3 + A(A + a)][2a 2 + 2aA - 3]\ 
^""iTTfA V 1_ V 3[2a + A] 2 J' 

In what follows we will analyze the stability of these critical points: 

• For point A: 

This fixed point is physically meaningless since = — ^f- < 0. 

• For point B, C: 

The existence of both points are only dependent on A. This can be understood as a 
result of the fact that the a-dependent term in Eq. (I17al) vanishes when — > 1 and 
z — > 0. Since the sign of \i\ is always opposite to the sign of /i 2 , these two points are 
saddle points. However, these a-independent critical points are found to be always stable 
and correspond to the future singularity in the standard FRW cosmology (30]. Therefore 
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the future singularity appearing in the standarad cosmology can be avoided by the loop 
quantum gravity effect. 

• For point D: 

This fixed point is physically meaningless since = — < 0. 

• For point E: 

The same as in the standard cosmology, this point is also unstable in LQC 

From the above analysis, we conclude that there is not any stable node for this kind 
of interacting dark energy model in LQC, and the future singularity appearing in the 



standard cosmology can be avoided. Our result o: 
of the numerical calculations obtained in Ref. 31] . 



analytical discussions agrees with that 



B. Case II: T = 3/?iJ(/^ + p m ) 



This type of interaction is motivated by analogy with dissipation of cosmological fluids 
and has been proposed for a possible dynamical solution to the coincidence problem [8|. 
The dynamics of this interacting phantom scalar field model in the standard FRW cos- 



mology has been studied in Ref. 



30j and it was found that there are two kinds of late 



time attractors. If the phantom field initially rolls down the potential the universe will ac- 
celerate forever and the total cosmic energy density decreases with the cosmic expansion, 
while if the phantom field initially climbs up its potential the universe will end with a big 
rip. Here we discuss the dynamics of this interacting model in LQC. For convenience, we 
introduce another variable £ = ^/Pm/ (V3H). The autonomous equations can be rewritten 
as: 

~ " 3 



x 

y' 

z' 



-3x 



V ^ 

2 y 2x(l-z) 



— x 



3x 2 - 



1 



-\l-Xxy-y 



3x z 



-3z - 3z(l - z){-x 2 
"1 ( , 1 



3 

2 VI 

y 2 



l 



i 2 2 

+ x — y 



1-z 
1 - 2z 



i 2 2 

+ x — y 



-3£ 



2M* 2 -2^ 2 



'1 - 2z) - 



2^(1 - z) 



;i - 2z) , (21a) 

(21b) 
(21c) 
(21d) 
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From Eq. (|T4j) it is easy to see that there are only three independent equations in the 
above system. This system has four critical points: 



Point A 



(22) 



xa 2 = ^(v / 1+ 4/3-1), y A = 0, 
£a 2 = ^1 + 4/3 + 1), z A = 0. 

Obviously this solution is physically meaningless since fi^ < 0. The solutions for the 
other three critical points are tedious and we do not present the details here. However, 
these solutions satisfy the following set of equations: 



P = f(x) , 



-x 2 - 



X 



+ 1 





2x z + 



(23a) 
(23b) 
(23c) 
(23d) 



where we have defined a cubic function 

f(x) = xi [ 2x + ^-j y-~^~ 

There is a critical point B with xb < 0, if 



(24) 



0<P<f 



-A - VA 2 + 12 
2V6 



(25) 



and this point B corresponds to an initially climbing- up phantom field [30J]. There are 
two other critical points [C and D) with xq,d > 0. One is physically meaningless and we 
label it by D. Point C exists for 



< (3 < min 



(26) 



and it corresponds to an initially rolling-down phantom field. 

Employing the standard techniques in the linearization and stability analysis, we can 
obtain three independent evolution equations of the linear perturbations. Using the eigen- 
values of the coefficient matrix for points B and C, found according to Eq. (|23|) . we find 
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that the critical point B is a saddle point, while point C is a stable critical point, and 
it is, therefore, a late time attractor. This is different from what was obtained in the 
standard FRW cosmology where both points B and C are late time attractors 30j. In the 
standard FRW cosmology, if the field initially rolls down the potential, the universe will 
enter a final state without Big Rip described point C, while for an initially climbing-up 
phantom field, the universe will enter a final state described by point B and end with a 
Big Rip. Therefore, the future singularity can be avoided by the loop quantum effect in 
LQC for the case in which the field initially climbs up the potential. As a result, in any 
case, there is no Big Rip in LQC. In Fig. (1) we show the stability regions for parameter 
space (A, 0). One can see that the effect of loop quantum gravity breaks the stability 
of the initially climbing-up field but leave that of the initially rolling-down field intact. 
This shows that if the universe does not evolve to a big rip in the standard cosmology its 
evolution seems to be uninfluenced by the loop quantum effect. 

Now, we will show results of numerical analysis we have carried out on this interacting 
phantom scalar field dark energy model in LQC. In Figs E]) we show the evolution 
curves of H and p for an initially climbing-up phantom field with different values of the 
coupling constant between the dark energy and dark matter. From these two figures 
we can see that H climbs up to reach its maximum value when the total cosmic energy 
density p reaches the value y = 0.75, and H goes down to zero when p reaches the 
maximum p = p c . After that, the universe contracts and then bounces. As time goes 
on, the universe will undergo oscillations with increasing frequency which may eventually 
blow up. Consequently, this seems to give rise to a new singularity. This kind of behavior 
of possible infinite frequency of oscillation also appears in the case of the interaction of 
the form V = a<j)p m , as found in Ref. {31]. In addition, the larger the value of /3, the later 
our universe enters the oscillating regime. 

In Figs. (|H El W) , we give the numerical results for an initially rolling-down phantom 
field with the requirement given in Eq. (l26j) . Fig. (j4|) shows the evolutionary properties 
of the universe with different initial conditions. Apparently the trajectories converge to 
the same final state determined only by parameter A. Figs. (0 E]) show the evolutionary 
curves of H and the effective equation of state for total cosmic energy w e ff. We find that 
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H does not oscillate and w e ff approaches finally to a constant which is less than — ^ but 
larger than —1. Thus the universe will keep accelerating forever while its energy density 
decreases with the cosmic expansion. Fig. ([7]) shows that in the final state the energy ratio 
of the phantom and dark matter reaches a constant. Therefore the coincidence problem 
can be alleviated. 

IV. CONCLUSION 

In conclusion, we have studied in loop quantum cosmology the dynamical system of a 
phantom field coupled to dark matter through an interaction of the form a<j)p m (case I) 
or 3/5iJ(p0 + p m ) (case II). The exponential potential for the phantom is used. For case I, 
there is a late time attractor solution in the standard FRW cosmology which corresponds 
to a big rip of the universe; whereas in LQC this solution transforms to be unstable, thus 
the big rip singularity which appears in the standard cosmology can be avoided by the 
loop quantum effect. 

For case II, it was found, in the standard FRW cosmology, that if the phantom field is 
initially rolling down the potential, the dynamical system has a late time attractor and 
the universe will accelerate forever, while if the field initially climbs up the potential, there 
is also a late time attractor but the universe will end with a big rip. By studying the 
dynamics of this interacting phantom model in LQC, we find that the universe with an 
initially climbing-up phantom field will oscillate forever; therefore the future singularity 
can be avoided by loop quantum effect; while for an initially rolling-down phantom field 
the universe will have the same late time evolution as that in the standard cosmology, 
i.e., the universe will accelerate forever with a constant ratio between the energies of the 
phantom and dark matter and the total energy density of the universe will decrease with 
the cosmic expansion. Therefore, in case II, the loop quantum effect only intervenes when 
the universe will evolve to a future singularity, and it sits idle when otherwise. 
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FIG. 1: The stability regions of (A, (3) parameter space for case II. In LQC, in the region I the 
rolling-down critical point (Point C) is a stable late time attractor; in region II there are no stable 
critical points. However, in the standard FRW cosmology, in the region I, both the climbing-up 
scaling solution and the rolling-down scaling solution are the stable late time attractors, while 
in the region II, the climbing-up solution is the stable late-time attractor. Ill represents the 
region of the solutions without physical meaning. 
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FIG. 2: The evolution of H for a climbing-up phantom field with different coupling constants 
{p = 0.2, 0.8) and A = 1, V = 1, p c = 1.5, O = -0.4 in case II. 




FIG. 3: The evolution of cosmic total energy density for a climbing-up phantom field with 
different coupling constants ((3 = 0.2,0.8) and A = 1, Vq = l,p c = 1-5, ^>o = —0.4 in case II. 
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FIG. 4: The convergence of different initial conditions to the attractor solution in the (x, y, z) 
phase space for a rolling down phantom field in case II. The red point denotes the critical point 
C 




FIG. 5: The evolution of H with time for a rolling-down phantom field with f3 = 0.2, A = 1,Vq = 
l,p c = 1-5 and 4>q = 0.4. 
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FIG. 6: The evolution of the effective equation of state for the total cosmic energy in a rolling 
down phantom model with (3 = 0.2, A = 1, Vq = 1, p c = 1.5 and 4>q = 0.4. 
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FIG. 7: The evolutionary curves of the fractional densities of a phantom field (dashed line) and 
dark matter (solid line) for a rolling-down phantom model with (3 = 0.2, A = 1, Vq = 1, p c = 1.5 
and 0o = 0.4. 
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